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ABSTRACT: High energy jet distributions measured since 1992 at the Fermilab Tevatron 
proton-antiproton collider are presented and compared to theoretical predictions. The statistical 
uncertainties on these measurements are significantly reduced relative to previous results. The 
systematic uncertainties are comparable in size to the uncertainty in the theoretical predictions. 
Although some discrepancies between theory and measurements are noted, the inclusive jet 
and dijet cross sections can be described by quantum chromodynamics. Prospects for reducing 
the uncertainty in the theoretical predictions by incorporating Tevatron measurements into 
the proton parton distributions are discussed. Dijet distributions, in excellent agreement with 
quantum chromodynamics, set a 2.5 TeV limit on the mass of quark constituents. 
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1 INTRODUCTION 

In the past seven years, quantum chromodynamics (QCD), the accepted theory of 
quark and gluon interactions, has been confronted with a set of precise and varied 
measurements of jet production from the Fermi National Accelerator Laboratory 
proton-antiproton collider. In pp collisions, jet production can be understood as 
a point-like collision of a quark or gluon from the proton and a quark or gluon 
from the antiproton. After colliding, these scattered partons manifest themselves 
as sprays of particles or "jets" . The extremely high energy of these interactions 
provides an excellent opportunity to test our understanding of perturbative QCD 
(pQCD). In particular, differences between QCD calculations and measured jet 
processes reveal information about both the partonic content of the proton and 
the nature of the strong interaction. Further, unexpected deviations may signal 
the existence of new particles or interactions down to distance scales of 10~^^ 
meters or less. 

Historically, jet measurements have involved tests of QCD at the highest ener- 
gies and searches for new physics. Prior to 1992, the inclusive jet cross section, 
one of the most important hadron collider results, had been measured at the ISR 
pp collider at a center-of-mass energy of ^/s = 63 GeV |^] , the CERN pp collider 
at 546 and 630 GeV [Q, |^, and the Fermilab Tevatron pp collider at 1800 GeV 
IQ, 1^. These data span a factor of twenty in beam energy and a factor of two 
hundred in jet energy transverse to the proton beam (Et) and, in general, are 
reasonably well described by QCD. 

Studies of the correlations between the leading two jets of an interaction also 
test QCD and provide opportunities to search for new physics. Measurements of 
dijet distributions prior to 1992, such as the invariant mass of the leading two 
jets or the angular separation between the jets, are also in good agreement with 
QCD §, 0, I, I, l^. Using the model of Eichten, Lane, and Peskin these 



data are used to search for composite quarks |§, ^ |5|, |^, |^, with the best 
limits coming from an inclusive jet cross section measurement |^. In this model 
excess jet production at large Et relative to QCD is interpreted as the product 
of quark constituent scattering. Because these data sets are relatively small all 
searches are limited by statistical precision rather than systematic uncertainties 
at the highest Et- A summary of these inclusive jet and dijet results can be 
found in the review of QCD by Huth and Mangano [T^ ]. 

In 1992, with the advent of high luminosity data collection at the Tevatron, 
a new era of precise pp jet physics began. From 1992-1993 (Run lA) each of 
the two collider experiments accumulated « 20 pb~^ of data and from 1994-1996 
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(Run IB) each accumulated ~ 100 pb~^. Both data sets were taken at ^/s = 1800 
GeV. At the end of Run IB a small portion of data ( ~ 600 nb~^) was taken at 
^/s = 630 GeV. The two runs together represent a data set ~ 50 times larger 
than any preceding data set. 

In contrast to previous tests of QCD, the Run 1 jet measurements have statisti- 
cal uncertainties significantly smaller than experimental or theoretical systematic 
uncertainties. In fact, the inclusive jet analysis from the Run lA sample yielded 
a significant discrepancy between data and contemporaneous theoretical predic- 
tions as well as good agreement with previous measurements [|l^. In addition 
to stimulating great excitement over the possibility of departures from QCD, the 
result motivated a re-evaluation of the uncertainties associated with inclusive jet 



cross section calculations |lj, 15, 1£]. Later measurements of the inclusive jet 



cross section using the Run IB samples stimulated further discussions since one 
confirmed [17| the Run lA measurement while the other was well described by 
QCD. In addition, subsequent studies have indicated that QCD can describe the 
observed Run lA cross section through adjustments to the parton distribution 
functions (PDFs), which represent the fraction of proton momentum carried by 
the constituent quarks and gluons. 

Independent of the actual high Et behavior of the inclusive jet cross section, 
the Run lA and IB results have revealed large uncertainties in the theoretical 
predictions. In particular, the PDFs are derived from fits to data from many 
different experiments, all of which are collected at low energy and extrapolated 
to Tevatron jet energies. Representation of the uncertainties in the resulting 
PDFs is a complex issue and has never been precisely resolved. Consequently, 
the Tevatron collaborations have taken an interest in derivation of the PDFs from 
information collected at the Tevatron. Measurements of dijet mass distributions 
and dijet cross sections at a variety of scattering angles provide information on 
the momentum distributions of the partons. By comparing the different mea- 
surements, constraints on the PDFs can be established, and sensitivity to the 
presence of new physics should ultimately be improved. 

This review opens with a description of the theoretical framework behind per- 
turbative QCD and a discussion of the uncertainties associated with predictions. 
Descriptions of the two detectors in which these measurements were performed 
and jet reconstruction algorithms are also provided. The review then devotes a 
section each to the Run 1 measurements of the inclusive jet cross section, the 
ratio of inclusive cross sections at different beam energies, dijet differential cross 
sections, and dijet mass and angular distributions. As will be seen, these mea- 
surements offer a detailed look into the composition of the proton and the nature 
of the strong interaction. A conclusion summarizes the results and suggests future 
avenues of research. 



2 PERTURBATIVE QUANTUM CHROMODYNAMICS 

2.1 Leading Order and N ext-to-Leading Order QCD 

The proton-antiproton interaction, a fairly general scattering process, nicely in- 
troduces the concepts of leading order and next-to-leading order pQCD. As 
illustrated in Fig. ||, inelastic scattering between a proton and an antiproton can 
be described as an elastic collision between a single proton constituent and single 
antiproton constituent. These constituents are collectively referred to as partons 
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and in QCD are quarks and gluons. The non-colliding constituents of the incom- 
ing proton and anti-proton are called beam fragments or spectators. Predictions 
for jet production are given by folding experimentally determined parton dis- 
tribution functions / with perturbatively calculated two-body scattering cross 
sections a . See, for example, reference |jl9| for a detailed discussion. The two 
ingredients can be formally combined to calculate any cross section of interest: 

o" = J2ij J dxidx2fi{xi, fi%)fj{x2, fiF)a[xiP, X2P, asinji), Q'^/i^f, Q'^ 1 1^^^- 

The parton distribution functions fk{x^[Xp) describe the momentum fraction x 
of the incident hadron momentum P carried by a parton of type k (gluons or 
quarks). The PDF is defined in terms of the factorization scale [xp. The hard 
two-body cross section is a function of the momentum carried by each of the 
incident partons xP, the strong coupling parameter a^, the scale Q characterizing 
the energy of the hard interaction, and the renormalization scale /x/j. Final state 
partons manifest themselves as collimated streams or "jets" of particles. This 
formal description includes no explicit hadronization or fragmentation functions 
to describe the transition from partons to jets. For most high Et measurements 
these effects are small and jets are identified as partons. 

The two-body scattering has been illustrated in Fig. || with a leading order 
(LO) graph. The cross section for this process is proportional to two powers of 
the strong coupling parameter, a^, which come from the two vertices. Although 
useful, the LO picture is too simple and has large normalization uncertainties. 
Next-to-leading order (NLO) or O(a^) calculations include one additional parton 
emission. A schematic example is shown in the second illustration of Fig. |^. Here 
a final state quark has radiated an additional gluon and the entire scattering 
process is proportional to . Depending on the proximity of the other partons, 
a "jet" could result from one or two (combined) partons. This in turn results 
in parton level predictions for the shape of jets, and for the effects of clustering 
parameters. 

A complete theoretical prediction of jet production should not depend on inter- 
nal calculational details; however, this is not the case for fixed order perturbative 
QCD, which depends on the choice of factorization and renormalization scales. 
The factorization scale, a free parameter, determines how the contributions of 
initial state radiation are factorized between the PDFs and a. The renormaliza- 
tion scale is also arbitrary and is related to choices in the theoretical calculations 
designed to control or renormalize ultraviolet singularities. Typically, as in this 
report, [iji and \if are set equal to each other, ^/j = [xp = /x, and /i is chosen 
to be of the same order as Q. The sensitivity of the theoretical predictions to \x 
is often taken as a measure of the uncertainty from the contributions of higher 
order terms. 

Predictions for jet production at NLO have been derived by Ellis, Kunszt, 
and Soper [^]; Aversa et al. |2l|; and Giele, Glover, and Kosower |2^. The 
NLO predictions are much less sensitive to the choice of /x [^^. The "EKS" 
program by Ellis, Kunszt, and Soper can generate analytic predictions for jet 
cross sections as a function of final state parameters. The "JETRAD" program 
by Giele, Glover, and Kosower |24] generates weighted "events" with final state 
partons. Cross sections are calculated by generating a large number of events as 
a function of final state parameters. All predictions in this document have been 
generated either with EKS or JETRAD. The two programs agree within a few 
percent [16 1 . 
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Figure 1: Factorization of the scattering process: In this example incoming quarks 
with momentum fraction xi and X2 of the incident hadrons scatter through gluon 
exchange and fragment into final state jets. 



2.2 Theoretical Choices 

A NLO QCD calculation requires selection of several input parameters including 
specification of a parton clustering algorithm, a perturbative scale, and PDFs. 
Taken together these choices can result in approximately 30% variations in the 
theoretical predictions. The single largest uncertainty is due to the PDFs. These 
input choices and the resultant changes in the predictions are discussed in the 
following paragraphs. The variations are illustrated by comparing predictions for 
the inclusive jet cross section. As in the experimental results to follow, the cross 
section is reported as a function of jet Et and pseudorapidity rj = — ln(tan(0/2)) 
where 9 is the angle between the jet and the proton beam, and is summed over 
(f), the azimuthal angle around the beam. For reference, rj = for jet production 
at 90 degrees relative to the proton beam. 

2.2.1 PARTON CLUSTERING 

The NLO predictions may include two or three final state partons. To convert the 
partons to the equivalent of jets measured in a detector, a clustering algorithm is 
employed. The Snowmass cone algorithm |2^ was proposed as a way to minimize 
the difference between theoretical predictions and experimental measurements. 
In theory calculations, two partons which fall within a cone of radius R in rj"(p 
space (R = \/ Ar]"^ + A(j)'^ and Arj and A(p are the separation of the partons in 
pseudorapidity and azimuthal angle) are combined into a "jet" . A standard cone 
radius of R = 0.7 is used by both collider experiments for most measurements. 
A consequence of this algorithm is that partons have to be at least a distance of 
2R apart to be considered as separate jets. 

Subsequent studies indicated that the experimental clustering algorithms (de- 
scribed in later sections) were more efficient at separating nearby jets p6| , p^ ] 
than the idealized Snowmass algorithm. In other words, two jets would be iden- 
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Figure 2: NLO theoretical predictions for the inclusive jet cross section. Each 
prediction is normalized relative to the NLO EKS calculation with CTEQ4M, 
/i = Et^^'/2, Rsep = 1.3. 

tified even though they were separated by less than 2R. An additional parameter, 
ep, was introduced in the QCD predictions to mimic the experimental effects of 
cluster in6rging and scpciration. Partons within Rsep 

X R were merged into a jet, 

while partons separated by more than Rsep 

X R were identified as two individual 
jets. As described in detail in the references, a value of Rsep = 1-3 was found 
to give the best agreement with the data for cross section and jet shape mea- 
surements [p^]. This corresponds in the data to the 50% efficiency point in jet 
separation p^], i.e. two jets within i?sep=l-3 of one another are merged 50% of 
the time and identified as two individual jets 50% of the time. At separations of 
l.OR the algorithms nearly always merge the two jets and at values of 1.6R they 
nearly always identify two separate jets. The top panel of Figure § shows the 
change in the prediction for the inclusive jet cross section using the EKS program 
when Rsep is decreased from 2 to 1.3. The result is primarily a normalization 
change of 5 to 7%. 

2.2.2 THE SCALE 

Because a NLO calculation is truncated at order there is some residual depen- 
dence on the scale // at which the calculation is performed. The scale is usually 
taken to be proportional to Et-'^^ (or just Et), maximum transverse jet energy 
^^max ^ given event, or the total center-of-mass energy. Other choices for 
the scale are possible. To study the scale dependence of the predictions, the 
magnitude of the scale is varied by a multiplicative coefficient, common choices 
are = 0.5Et^''^ and 0.5£;t The middle panel of Fig. | shows a mild Et 
dependence at the 2-9% level on the definition of the scale. The lowest panel 
shows that typical variations of the multiplicative coefficient lead to 5-20% shifts 
in the cross section with only small changes in shape above 100 GeV. 
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2.2.3 PARTON DISTRIBUTION FUNCTIONS 

Parton distribution functions are derived from global fits to data primarily from 
deep inelastic scattering (DIS) fixed target and electron-proton collider experi- 
ments. In DIS, a lepton probe is used to sample the partonic structure of the 
target hadrons. For use at the Tevatron, the resulting distributions must be 
evolved from the low energy DIS results to the high energy range of jet measure- 
ments. 

As the data from DIS experiments increases and improves, new parton distri- 
bution functions are derived. Some recent PDFs also incorporate data from the 
Tevatron. The result is a plethora of PDFs, each with its own specific list of 
data included in the fit, assumptions as to the value of as, functional forms for 
the quark and gluon momentum distributions, and assumptions concerning the 
contributions of the gluons, which are not constrained by the DIS results. The 
most recent PDFs include the most precise data and combined knowledge, and 
supplant the previous PDFs. 

Several groups have analyzed the available data and produced families of can- 
didate PDFs. Within a family, the individual PDFs represent the range of pre- 
dictions resulting from changes in one of the input assumptions such as the value 
of as- Additional variations come from differences in the input data sets or the 
relative weighting between the data sets. The CTEQ2M [||, CTEQ3M [||] and 
MRSA' families of PDFs incorporate data published before 1994 and do not 
include Tevatron jet data. The CTEQ4M PDFs include data published be- 
fore 1996; CTEQ4HJ |^l|] additionally has a high x gluon adjustment designed to 
accommodate the Run lA high Et jet cross section measurements. The MRST 
1 32 1 PDF family utilizes data published before 1998 and adds a contribution for 



a putative initial transverse momentum of the partons. 

Comparison between observed cross sections and NLO predictions with alter- 
nate PDFs provide some insight into the quark and gluon composition of the 
proton. As shown in the Fig. ^ the PDFs can result in ~ 20% variations in jet 
cross sections; typical variations within a set are on the order of 5-10%. Although 
the quark distributions were thought to be well known and the contribution of 



the gluons was expected to be small, investigations [33, 34 1 have shown that there 



are uncertainties at high x, which were ignored in the derivation of early PDFs. 



In addition, studies ||3l|] revealed that the gluon distribution could be adjusted to 
give a significant increase in the jet cross section at high Ej-, while maintaining 
reasonable agreement with the low energy data sets. 



3 DETECTORS AND JET DEFINITIONS 



3.1 The CDF and DO Detectors 

The Collider Detector at Fermilab (CDF) [H and the DO Detector Hi at the 
Fermi National Accelerator Laboratory Tevatron pp collider are complementary 
general purpose detectors designed to study a broad range of particle physics top- 
ics. The Tevatron is designed with six locations at which the counter-rotating 
proton and antiproton beams can collide. The CDF and DO detectors each oc- 
cupy one of these interaction regions; DO takes its name from the alphanumeric 
designation of its interaction region. Each detector is comprised of a series of 
concentric sub-systems, which surround the pp interaction region. Immediately 
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Figure 3: NLO theoretical predictions for the inclusive jet cross section using dif- 
ferent PDFs. Each prediction is normalized relative to the NLO EKS calculation 
with CTEQ4M. All used i_i = Et/2 and Rsep = 1.3 



outside the Tevatron beam pipe each detector has a tracking system designed to 
detect charged particles. The tracking systems are surrounded by calorimeters, 
which measure the energy and direction of electrons, photons, and jets. Finally, 
the calorimeters are covered with muon tracking systems. 

Luminosity monitors at each detector measure the beam exposure. These mon- 
itors are located near the beam lines and detect particles from the pp collisions. 
The luminosity at a given site is given by dividing the luminosity monitor event 
rate by that portion of the total pp cross section to which the luminosity monitors 
have acceptance. A direct comparison of jet cross sections at the two experiments 
is complicated by the fact that the collaborations adopt slightly different values 
for the total pp cross section [37, 39 1. As a result, the CDF jet cross sec- 
tion measurements are 2.7% higher than the corresponding DO jet cross section 
measurements |37]. 

In the following sections we mention primarily the central tracking systems and 
the calorimeters since these systems are the most important for jet identification 
and reconstruction. As will be seen, CDF has a high resolution particle tracking 
system, which makes a crucial contribution to jet energy calibration. On the 
other hand, and in a complementary fashion, DO has highly segmented, uniform, 
and thick calorimetry well suited to in-situ jet energy measurements. 

Both calorimeters are segmented into projective towers. Each tower points 
back to the center of the nominal interaction region and is identified by its pseu- 
dorapidity and azimuth. The polar angle 9 in spherical coordinates is measured 
from the proton beam axis, and the azimuthal angle (p from the plane of the 
Tevatron. The towers are further segmented longitudinally into individual read- 
out cells. The energy of a calorimeter tower is obtained by summing the energy 
in all the cells of the same pseudorapidity and azimuth. The transverse energy 
or Et of a tower is the sum of the Et components of the cells. These compo- 
nents are calculated by assigning each cell of a tower a massless four-vector with 
magnitude equal to the energy deposited and with the direction defined by the 
unit vector pointing from the event origin to the center of the tower segment. 
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Figure 4: The distribution of energy in a dijet event in the DO detector. 
3.1.1 DO in Run 1 

The central tracking volume (|r/| ~ 2) of DO includes the vertex chamber, tran- 
sition radiation detector, and central drift chamber, arranged in three cylinders 
concentric with the beamline, and two forward drift chambers. The non-magnetic 
tracker is compact, with an outer radius of 75 cm and an overall length of 270 cm 
centered on z = 0. Without the need to measure momenta of charged particles, 
the prime considerations for tracking were good two track resolving power, high 
efficiency, and good ionization energy measurement. For jet physics the central 
tracker is used to find event vertices. 

Jet detection in the DO detector primarily utilizes liquid argon-uranium 
calorimeters which are hermetic, finely segmented, thick, uniform and have unit 
gain. These sampling calorimeters are composed of alternating layers of liquid 
argon and absorber. The particles in a jet interact with the absorber, and the 
resulting particle shower ionizes the liquid argon. This ionization is detected as 
a measure of the jet energy. The calorimeter is enclosed in three cryostats: the 
central calorimeter (CC) covers \7]\ < 1.2, and the end calorimeters (ECs) extend 
the coverage to \ri\ < 4.4. The calorimeters have complete 27r azimuthal cover- 
age. Between the cryostats the calorimeter sampling is augmented by scintillator 
tiles with segmentation matching the argon calorimeters. Fig. ^ is a graphical 
representation of the calorimeter modules described below. 

The CC includes three concentric rings of calorimeter modules. There are 32 
electromagnetic calorimeter modules (CCEM) in the inner ring, 16 fine hadronic 
modules (CCFH) in the middle ring, and 16 coarse hadronic modules (CCCH) in 
the outer ring. The CCEM and CCFH calorimeters have uranium absorber plates, 
and the CCCH has copper absorber plates. Longitudinal segmentation includes 
seven samples: four in the CCEM, three in the CCFH, and one in the CCCH. At 
rj = the CC has a thickness of 7.2 nuclear absorption lengths. The calorimeter 
cells are segmented into Arj x A(p = 0.1 x 0.1 except at shower maximum in the 
third layer of the CCEM where the segmentation is 0.05 x 0.05. The calorimeter 
segmentation is designed to form projective towers of size Ar] x Acp = 0.1 x 0.1 
geometry which point to z = 0. In the CC because of the finer resolution in the 
third layer of the CCEM each tower is comprised of eleven cells. 
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The two mirror-image ECs contain four module types. An electromagnetic 
module (ECEM) surrounding the beam pipe is backed by an inner hadronic mod- 
ule (ECIH) which also surrounds the pipe. Outside the ECEM and ECIH are 
concentric rings of 16 middle and outer hadronic modules (ECMH and ECOH). 
The ECEM and ECIH have uranium absorber plates and four longitudinal seg- 
ments. The ECMH has four uranium absorber segments nearest the interaction 
point and one stainless steel absorber segment behind. The ECOH employs 
stainless steel plates and includes three longitudinal segments. At |?7|= 2.0 a 
jet encounters 8-9 longitudinal segments and a nuclear interaction thickness of 
~ IOAq. The transverse segmentation is similar to that of the CC. The CC(EC) 
electron energy resolution is 14.8(15.7)%/\/ E plus 0.3(0.3)% added in quadrature 
and for hadrons 47.0(44.6)%/\/!E plus 4.5(3.9)% added in quadrature. 

The calorimeter response to the different types of particles is the most im- 
portant aspect of the DO jet energy calibration. Electromagnetically interact- 
ing particles, like photons (7) and electrons deposit most of their energy in the 
electromagnetic sections of the calorimeters. Hadrons, by contrast, lose energy 
primarily through nuclear interactions and extend over the full interaction length 
of the calorimeter. In general, the calorimeter response to the electromagnetic 
(e) and nuclear or hadronic components (/i) of hadron showers is not the same. 
Non-compensating calorimeters have a response ratio e/h greater than one and 
suffer from non-Gaussian event-to-event fluctuations in the fraction of energy 
lost through electromagnetic production. Such calorimeters give a non-Gaussian 
signal distribution for hadrons and jets. The DO calorimeter is nearly compen- 
sating and the hadronic and jet response are well described by a Gaussian dis- 
tribution pO[ |. In fact, single jet resolution as measured with dijet events has a 
Gaussian line-shape and is approximately 7% at 100 GeV and 5% at 300 GeV 

0- 

The DO jet calibration requires correction for the hadronic response of the jet; 
showering of energy outside the cone; and subtraction of an offset which can 
be attributed to instrumental effects, pile-up from previous beam-beam cross- 
ings, additional interactions, and spectator energy. These corrections are derived 
primarily in-situ [^]. The correction for hadronic response begins with the elec- 
tromagnetic calibration of the calorimeter which is performed with dielectron 
and diphoton decays of the Z and vr" resonances since electrons and photons de- 
posit their energy in the electromagnetic calorimeters. The hadronic response for 
centrally produced 7-jet events is derived from data using Et balance between 
the photon and jet. Because of the rapidly falling 7-jet cross section the central 
calorimeter balance technique is limited to jet energies below 150 GeV. Balancing 
central, well measured photons against high rapidity jets permits the energy cali- 
bration to exceed 300 GeV. The response calibration is extended above 300 GeV 
using simulated 7-jets events. For a R=0.7, 100 GeV jet at rj = the hadronic 
response is about 0.85 it 0.01. The showering correction for this jet is about 
1 lb 1% and is derived from a study of jet profiles. The total offset correction 
for a R=0.7 jet at r/ = as determined from a study of data sets with various 
triggering and luminosity conditions is about 2.6 it 0.3 GeV. The contribution 
from the underlying event or spectator energy alone is roughly 0.9 it 0.1 GeV. At 
r] = the mean total energy correction for a R=0.7, 100 GeV jet is 15.0 it 1.7% 
and decreases slowly with energy. 
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3.1.2 CDF in Run 1 

The CDF central (|r/| ~ 1) tracking system has three sub-systems located within 
a 1.5 T magnetic field, that is provided by a superconducting solenoid coax- 
ial with the beam. Nearest the beam, a four-layer silicon microstrip vertex 
(SVX) detector occupies the radial region between 3.0 and 7.9 cm from 
the beamline and provides precision r — (j) measurements. Outside the SVX a 
vertex drift chamber (VTX) provides r-z tracking information and is used to 
locate the position of the pp interaction (event vertex) in z, along the beam line. 
Both the SVX and the VTX are mounted inside a 3.2 m long drift chamber 
called the central tracking chamber (CTC). The radial coverage of the CTC is 



from 31 to 132 cm. The momentum resolution [35| of the SVX-CTC system is 
5Pt/P^ = [(0.0009Pt)^ + (0.0066)2]V2 ^here Pt has units of GeV/c. The CTC 
provides in-situ measurement of the calibration and response of the calorimeter 
to low energy particles (where test beam information is not available) along with 
measurements of jet fragmentation properties. 

The solenoid and tracking volumes of CDF are surrounded by calorimeters, 
which cover 27r in azimuth and |r/| < 4.2. The central electromagnetic (CEM) 
calorimeter covers \r]\ < 1.1 and is followed at larger radius by the central hadronic 
calorimeters (CHA and WHA) which cover \t]\ < 1.3. These calorimeters use 
scintillator as the active medium. The CEM absorber is lead and the CHA/WHA 
absorber is iron. The calorimeters are segmented into units of 15 degrees in 
azimuth and 0.1 pseudorapidity. Two phototubes bracket each tower in </> and 
the average of the energy in the two tubes is used to determine the (j) position 
of energy deposited in a tower. The interaction length of both the CHA and 
WHA is 4.5Ao. Electron energy resolution in the CEM is 13.7%/VE plus 2% 
added in quadrature. For hadrons the single particle resolution depends on angle 
and varies from roughly 50% /^/E plus 3% added in quadrature in the CHA 
to 75%) /VE plus 4% added in quadrature in the WHA. In the forward regions 
(1.1 < \r]\ < 4.2) calorimetric coverage is provided by gas proportional chambers. 
The plug electromagnetic (PEM) and hadronic calorimeters (PHA) cover the 
region from 1.1 < \rj\ < 2.4. The forward electromagnetic (FEM) and hadronic 
calorimeters (FHA) cover the region from 2.2 < |?7| < 4.2. The segmentation of 
these detectors is roughly 0.1 in r] and 5 degrees in (j). 

Figure ^ shows a 3-jet event in the CDF calorimeter. In this "lego" plot the 
calorimeter is "rolled out" onto the r]-(f) plane. The tower height is proportional 
to the Ej- deposited in the tower. The darker and lighter shading of each tower 
corresponds to the Ej- of the electromagnetic and hadronic cells of the tower 
respectively. The oval around each clump of energy indicates the jet clustering 
cone described in the next section. 

Observed jet energies are corrected for a number of effects including the calibra- 
tion and response of the detectors and the background energy from the remnants 
of the pp interaction. In the central detectors corrections associated with de- 
tector response are obtained from a Monte Carlo program, which is tuned to 
give good agreement with the data from electron and hadron test beams {Ebeam 
ranged from 10 to 227 GeV), and to data from an in-situ study using CTC track 
momenta and isolated hadrons with 400 MeV/c< Pt < 10 GeV/c. In addition, 
the Monte Carlo jet fragmentation parameters were tuned to agree with both 
the number and momenta of particles observed in jets. The resulting simulation 
was then used to determine "response functions", which represent the response 
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of the calorimeter to jets as a function of jet Et- The width of the response 
functions represents the jet energy resolution and can be expressed as a = O.IEt 
+ 1 GeV m. Calibration of the plug and forward detectors is achieved through a 
dijet balancing technique |43]. The of jets in the plug and forward detectors 
is balanced against jets in the central region whose calibration is pinned by the 
tracking information. The energy left in a jet cone from the hard interaction spec- 
tators is called the "underlying event energy" . A precise theoretical definition of 
this quantity does not exist. Experimentally, it has been estimated from events 
in which there is no hard scattering, and from the energy in cones perpendicular 
to the jet axis in dijet events. For a cone of radius 0.7 the contribution to the 
jet Et is of order 1 GeV. An uncertainty of 30% is assigned to this quantity to 
cover reasonable variations in the definition of this quantity. The magnitude of 
the total jet energy corrections and the corresponding uncertainty depend on the 
Et of the jets and their location in the detector, and on the specific analysis. 
Typically, the jet Et correction factor is in the range of 1.0 to 1.2. 



3.1.3 CDF and DO in Run 2 

The Fermilab Tevatron is undergoing major improvements for Run 2, which is 
slated to begin in the year 2000. The energy of the proton and antiproton beams 
will increase from 900 to 1000 GeV and the instantaneous luminosity will increase 
by at least a factor of 5. The expected Run 2 data sample is 2fb~^ in the first 2 
years of operation. 

Both the CDF and DO collaborations are improving their detectors for Run 
2 operations. Among other things, CDF is replacing the gas calorimeters with 
scintillating tiles and is closing gaps between calorimeters. CDF is also extending 
the tracking such that calibration of the new plug detector can be accomplished 
using tracks, as in the central detector. DO is also upgrading a number of detector 
components. For jet measurements the most important upgrade is the replace- 
ment of the nonmagnetic tracking system with a magnetic tracker. The tracker 
(including a high resolution silicon strip detector, a scintillating fiber tracker, and 
electron/photon preshower detectors) will enable calibration of the calorimeter 
with single particles. Thus, in Run 2 the detectors will be more similar than 
in Run 1: CDF will have more uniform calorimetry and DO will have magnetic 
tracking. The result of these improvements should be an overall reduction of jet 
measurement uncertainties. 



3.2 The Experimental and Theoretical Definitions of a Jet 

From the experimental point of view, jets are the manifestation of partons as 
showers of electromagnetic and hadronic matter. Jets are observed as clusters 
of energy located in adjacent detector towers. Typically, a jet contains tens of 
neutral and charged pions (and to a lesser extent kaons), each of which showers 
into multiple cells. A single jet illuminates roughly 20 towers (this corresponds 
to approximately 40 calorimeter cells in the CDF detector and over one hundred 
detector cells in the DO detector) . Fig. ^ illustrates the energy distribution of a 
two-jet event in the DO detector. Each rectangular outline represents the energy 
deposited in calorimeter cells at fixed t] and depth. The eye quite naturally 
clusters the energy into two jet-like objects. However, for the purposes of jet 
cross section measurements a more quantitative definition of the jet is required. 
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For high Ej- measurements at CDF and DO, cone algorithms are used to 
identify jets in the calorimeters. Cone algorithms operate on objects in pseu- 
dorapidity and azimuth space, such as particles, partons, calorimeter cells, or 
calorimeter towers. For the sake of simplicity towers will be used in the following 
discussions. The CDF and DO algorithms are both based on the Snowmass algo- 
rithm pH]. This algorithm defines a jet as those towers within a cone of radius 



AR = \J A(/)^ + Ary^ = 0.7. The jet is the sum of the transverse energies of 
the towers in the cone, and the location of the jet is defined by the E^ weighted 
77 and ^ centroids: 

Ei^' = Y,ET' 

i 

Vjet = {J2ET\)/ET^'^' 



The sum over i is over all towers that are within the jet radius R. 

In Figure |5| the clustering cone is shown by the oval around each jet and 
is centered on the cluster centroid. The two overlapping cones in this event 
indicate that the two nearby clusters have been identified as two separate jets. 
The Snowmass algorithm doesn't specify cell thresholds or the handling of such 
overlapping jet cones. These details must be dealt with according to the needs 
of individual experiments, and led to the introduction of the Rsep parameter in 
the theoretical calculations. 



In the DO experiment, jets are defined in two stages |27]. In the first or cluster- 
ing stage, all the energy that belongs to a jet is accumulated, and in the second 
stage, the r], (j) and Et of the jet are defined. The clustering consists of the fol- 
lowing steps: 1) Calorimeter towers with Et > 1 GeV are enumerated. Starting 
with the highest Et tower, preclusters are formed by adding neighboring towers 
within a radius of R=0.3. 2) The jet direction is calculated for each precluster 
using the sums defined for the Snowmass algorithm. 3) All the energy in the 
towers in a cone of radius R=0.7 around each precluster is accumulated and used 
to recalculate r] and (j). 4) Steps 2 and 3 are repeated until the jet direction is 
stable. Finally, in the second stage, the jet energy and directions are calculated 
according to the equations: 



{i:^EiY + {Y.^EiY 

tan 6 jet = ^ ^. 

jet 

<Pjet = i^n-\Y.Ei/Y.K? 
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where E^. = Ei siniOi) cos(0i), Ey = Ei sin(0j) sin(0j), and E\ = Ei cos{6i). Stud- 
ies have shown that for rj <~ 1 the final DO jet directions are, within experimental 
errors, equal to the Snowmass directions. Overlapping jets are merged if more 
than 50% of the smaller jet E^ is contained in the overlap region. If less than 
50% of the energy is contained in the overlap region, the jets are split into two 
jets and the energy in the overlap region is assigned to the nearest jet. After 
merging or splitting the jet directions are recalculated. 



The CDF cluster algorithm [44| has two stages of jet identification similar to 
the DO algorithm. The first stage consists of the following steps: 1) a list of 
towers with Et >1.0 GeV is created; 2) preclusters are formed from an unbroken 
chain of contiguous seed towers with continuously decreasing tower Et', if a tower 
is outside a window of 7 x 7 towers surrounding the seed it is used to form a new 
precluster; 3) the preclusters are grown into clusters by finding the Ej- weighted 
centroid and collecting the energy from all towers with more than 100 MeV 
within R=0.7 of the centroid; 4) a new centroid is calculated from the set of 
towers within the cone and a new cone drawn about this position; steps 3 and 4 
are repeated until the set of towers contributing to the jet remains unchanged; 5) 
overlapping jets are merged if they share >75% of the smaller jet's energy; if they 
share less, the towers in the overlap region are assigned to the nearest jet. In the 
second stage, the final jet parameters are computed. The angles and energy are 
calculated as in the DO algorithm, however the jet Et is given by: 

Et = EsinO jet- 
One philosophical difference between the CDF and Snowmass algorithms is that 
the CDF jets have mass {Et 7^ Pt)- Studies |^ found that the CDF clustering 
algorithm and the Snowmass algorithm were numerically very similar. 

Both the CDF and DO algorithms are based on the Snowmass algorithm, but 
the reality of jets of hadrons measured by a finitely segmented calorimeter necessi- 
tates the introduction of additional steps and cuts. Apart from small definitional 
details, the most significant difference between the Snowmass algorithm and its 
experimental implementations is the handling of cluster merging and separation. 
To simulate these effects in the NLO calculations the Rsep parameter was in- 
troduced. Considering the complexity of the hadronic showers is it remarkable 
that a NLO calculation, with only 2 or 3 partons in the final state and a single 
additional parameter can provide a good description of the observed jet shapes 



and cross sections [26|. 



4 THE INCLUSIVE JET CROSS SECTION AT 1800 GEV 



4.1 Introduction 

The inclusive jet cross section represents one of the most basic tests of QCD 
at a hadron-hadron collider. It reflects the probability of observing a hadronic 
jet of a given Et and rapidity in a pp collision. The term inclusive indicates 
that all the jets in an event are included in the cross section measurement and 
that the presence of additional non-jet objects (for example electrons or muons) 
is irrelevant. Theoretical calculations are normally expressed in terms of the 
invariant cross section: 



Ed^a/dp^ 
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where E and p are the jet energy and momentum. In terms of the natural 
experimental variables, the cross section is given by 

d'^a/dErdr] 

which is related to the first expression by, 

Ed^a/dp^ = {l/2T:ET)d^a/dETdr]. 

Massless jets have been assumed {Pt = Et); the 2tt comes from the integration 
over the azimuthal angle (p. The measured cross section is simply the number 
of jets observed in an r/ and Et interval normalized by the total luminosity 
exposure, C: 

d^a/dErdri = N/{AEtAt]C). 

The CDF and DO inclusive jet analyses place similar requirements on the events 
and jets selected for calculation of the cross sections. Both collaborations elimi- 
nate poorly measured events by requiring the event vertices to be near the center 
of the detector. Backgrounds (such as cosmic rays) are eliminated by rejecting 
events with large missing Et- Spurious jets from cosmic ray or instrumental 
backgrounds are eliminated with quality cuts based on jet shapes. Corrections 
are made to the measured cross sections to account for the event and jet detection 
inefficiencies, mismeasurement of the jet energies, and for energy falling in a jet 
cone from other sources (such as remnants of the pp collision, or additional pp 
interactions). No corrections are made for partons showering outside the jet cone 
as this should be included in the NLO theoretical calculations. 

We begin our discussion of the inclusive jet cross section with measurements by 
the CDF collaboration. The Run lA measurement stimulated interest due to a 
discrepancy with QCD predictions at high Et- A subsequent measurement by the 
DO collaboration with Run IB data is well described by theoretical predictions, 
while the preliminary CDF Run IB result is in good agreement with the lA 
measurement. In the next sections, the CDF and DO measurements are each 
described in some detail and then comparisons of the two results are presented. 



4-2 The CDF Cross Section 

In 1996 CDF published the inclusive jet cross section measured from the Run 



lA data sample [13| for jet Et from 15 to 440 GeV in the central rapidity range 
0.1< \r]\ < 0.7. This analysis followed the same procedures as previous measure- 
ments ^ for correcting the cross section and for estimating the systematic 
uncertainties. This process will be described briefly below. More detail is avail- 
able in references @, 13, |45|| . 



The measured jet Et spectrum requires corrections for energy mismeasurement 
and for smearing effects caused by finite Et resolution. This is accomplished 
with the "unsmearing procedure" described in [45|. As mentioned earlier, a 



Monte Carlo simulation program [g] tuned to the CDF data is used to determine 
detector response functions. A trial true (unsmeared) spectrum is smeared using 
these response functions and compared to the raw data. The parameters of the 
trial spectrum are iterated to obtain the best match between the smeared trial 
spectrum and the raw data. The corresponding unsmeared curve is referred to 
as the "standard curve", and is used to correct the measured spectrum. The 
simultaneous correction for response and resolution produces a result which is 
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Figure 6: Inclusive jet cross section measured by CDF using the Run lA data 
sample compared to predictions from the EKS program with /i = £'jf*/2, Rsep 
= 2.0 and MRSDO'. The error bars represent the statistical uncertainty on the 
cross section. The shaded band represents the quadrature sum of the correlated 
systematic uncertainties. Also shown is the percentage difference between the 
predictions using other PDFs and the prediction with MRSDO'. 



independent of the Et binning used in the measurement while preserving the 
statistical uncertainty on the measured cross section. For jet Et between 50 and 
about 300 GeV these corrections increase both the Et and the cross section by 
~10%. At lower and higher Et the cross section corrections are larger due to the 
steepening of the spectrum. 

Figure ^ shows the CDF measurement of the inclusive jet cross section from 
the 20 pb~^ Run lA data sample compared to a NLO prediction from the EKS 
program with = Et/2, Rsep = 2.0 and MRSDO' PDFs. The inset shows the 
cross section on a logarithmic scale. There is good agreement between data and 
theory over eleven orders of magnitude. The main figure shows the percentage 
difference between the data (points) and theory. The bars on the points represent 
the statistical uncertainty. While excellent agreement is observed below 200 GeV, 
an excess of events over these theoretical predictions is observed at high Et- 
Predictions with other PDFs which were available at the time are shown by the 
additional curves. For these predictions the percentage from the default theory 
(MRSDO') is shown. As can be seen from the figure, the best agreement below 
200 GeV is with the MRSDO' PDF. Other PDFs agree less well, and none show 
the rise at high Et observed in the data. 

The systematic uncertainties on the Run lA cross section were evaluated using 
the procedures in reference p^]. In short, new parameter sets are derived for 
±1 standard deviation shifts in the unsmearing function for each source of sys- 
tematic uncertainty. The parameters for the Run lA cross section are given in 
Reference [13|. Figures ^(a-h) show, for the Run lA data sample, the percentage 
change from the standard curve as a function of Et for the seven largest system- 
atic uncertainties: (a) charged hadron response at high Pt; (b) the calorimeter 
response to low-P^ hadrons; (c) ±1% on the jet energy for the stability of the 
calibration of the calorimeter; (d) jet fragmentation functions used in the simula- 
tion; (e) ±30% on the underlying event energy in a jet cone; (f) detector response 
to electrons and photons and (g) modeling of the detector jet energy resolution. 
An eighth uncertainty, an overall normalization uncertainty of ±3.8%, was de- 



18 



Gerald C. Blazey and Brenna L. Flaugher 



Energy Scale Stability 1 ^ 



Underlying Event 30% 



Calorimeter Resolution 10% 



100 200 300 

Transverse Energy (CieVJ 



Low P-p Hadron response 5% 



Fragmentation 



Electron/Photon Response 2% 



Normalization 3.8% 



lOO 200 300 400 

Transverse Energy (CieV> 



Figure 7: The percentage change in the Run lA inclusive jet cross section when 
various sources of systematic uncertainty are changed by ±1 standard deviation 
from their nominal value. 



rived from the uncertainty in the luminosity measurement pq] (±3.5%) and the 
efficiency of the acceptance cuts (±1.5%). The eight uncertainties arise from 
different sources and are not correlated with each other. The la shifts are eval- 
uated by changing only one item at a time in the Monte Carlo simulation, such 
as high pt hadron response. The resulting uncertainty in the unsmeared cross 
section is thus 100% correlated from bin to bin, but independent of the other 
seven uncertainties. 

To analyze the significance of the Run lA result, the CDF collaboration used 



four normalization-independent, shape-dependent statistical tests [|13|. The eight 
sources of systematic uncertainty are treated individually to include the Et de- 
pendence of each uncertainty. The effect of finite binning and systematic uncer- 
tainties are modeled by a Monte Carlo calculation. Between 40 and 160 GeV, the 
agreement between data and theory is >80% for all four tests. Above 160 GeV, 
however, each of the four methods yields a probability of 1% that the excess is 
due to a fluctuation. If the test is performed for other PDFs, agreement at low 
Et is reduced, as is the significance of the excess at high Et- The best agreement 



at high Et for the curves shown is with CTEQ2M[28| which gives 8%, but the 
low Et agreement is reduced to 23%. The excess of events observed at high Et 
initiated a re-evaluation of the uncertainties in the PDFs, particularly at high x. 
One outcome of this re-examination was CTEQ4HJ. This PDF incorporates the 
same low energy DIS data as in CTEQ4M, but the high Et jet data is weighted 
to accentuate its contribution to the global parton distribution fit in the high Q 
and X region. 



Figure ^ shows the preliminary Run IB result |47| compared to the Run lA 
results and to CTEQ4M (a more modern PDF than the ones in Figure I). The 
Run lA and IB data sets are in excellent agreement. Only statistical uncertainties 
are shown. The greatly reduced statistical error bars on the Run IB data are 
due to the five-fold increase of luminosity between the two runs. Analysis of 
the Run IB data follows exactly the sequence of the Run lA data with some 
additional corrections specific to the Run IB running conditions. The Run IB 
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Figure 8: The preliminary Run IB inclusive jet cross section compared to the 
published Run lA data and to a QCD prediction from the EKS program with 
^= Rsep, = 1.3 and CTEQ4M. 



systematic uncertainties are very similar to the uncertainties derived for the Run 
lA sample. Figure ^ shows the Run IB result compared to predictions with other 
recent PDFs. Although CTEQ4HJ is a bit contrived the good agreement with 
the data demonstrates the flexibility of the PDFs and the ability of the QCD 
calculations to describe these Run lA and Run IB inclusive jet cross sections. 
Quantitative comparisons between the Run IB data and theoretical predictions 
are underway. 



4-3 The DO Cross Section 

In 1998 the DO collaboration finalized a 92 pb"^ measurement of the inclusive 



jet cross section |18]. The DO analysis differs from the CDF analysis in that 
the spectrum is corrected independently for energy calibration and then for dis- 
tortion by finite jet energy resolution. After passing the various jet and event 
selection criteria each jet is corrected individually for the average response of the 
calorimeter. As mentioned in Section 3 the jet response was determined with well 
measured photon-jet events. However, the background free, energy corrected Et 
spectrum still remains distorted by jet energy resolution. The distortion is cor- 
rected by first assuming that an ansatz function {AEt~^) ■ (1 — 2Et / \/s)'-' will 
describe the actual E^ spectrum, then smearing it with the measured resolution, 
and finally comparing the smeared result with the measured cross section. The 
procedure is repeated by varying parameters A,B, and C until the best fit is 
found between the observed cross section and the smeared trial spectrum. At all 
Et, the resolution (measured by balancing Ej- in jet events) is well described by 
a Gaussian distribution; at 100 GeV the standard deviation is 7 GeV. The ratio 
of the initial ansatz to the smeared ansatz is used to correct the cross section 
on a bin~by-bin basis |^]. The resolution correction reduces the observed cross 
section by (13±3)% [(8±2)%] at 60 GeV [400 GeV]. 

Fig. Iin shows the final inclusive jet cross section as measured by the DO col- 



laboration in the rapidity region |7y| < 0.5 [18|. This rapidity interval was chosen 
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Figure 9: The percent difference between the prehminary Run IB inclusive jet 
cross section and QCD predictions from the EKS program with fj,= Rsep 
= 1.3, and a variety of current PDFs. 



Table 1: DO Inclusive jet cross section total uncertainty correlations. 
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since the detector is uniformly thick (seven or more interaction lengths with no 
gaps) and both jet resolution and calibration are precise. The figure shows a theo- 
retical prediction for the cross section from JETRAD . There is good agreement 
over seven orders of magnitude. For the calculation shown here /i = 0.5£'j?^^, the 
PDF is CTEQ3M, and 7^sep=1.3. 

Fig. 11 shows the cross section uncertainties. Each curve represents the average 
of nearly symmetric upper and lower uncertainties. The energy scale uncertainty 
which varies from 8% at low Et to 22% at 400 GeV dominates all other sources 
of uncertainty, except at low Et, where the 6.1% luminosity uncertainty is of 
comparable magnitude. The total systematic error is 10% at 100 GeV and 23% 
at 400 GeV. Although the individual errors are independent of one another each 
error is either 100% or nearly 100% correlated point-to-point, and the overall 
systematic uncertainty is highly correlated. Table |l| shows that the bin-to-bin 
correlations in the full uncertainty for representative Et bins are greater than 
40% and positive. The high degree of correlation will prove a powerful constraint 
on data-theory comparisons. 

Figure shows the ratios {D-T)/T for the DO data (D) and JETRAD NLO 
theoretical predictions (T) based on the CTEQ3M, CTEQ4M and CTEQ4HJ for 
the region \r]\ < 0.5. Figure ^ shows the same ratios for the 0.1 < |?7| < 0.7 
data. Given the experimental and theoretical uncertainties, the predictions are 
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Table 2: comparisons between JETRAD and DO , |r/| < 0.5 and 0.1 < \rj\ < 
0.7 data for fj, = O.SE'™^^, 7^scp=l-37^, and various PDFs. There are 24 degrees 
of freedom. 



PDF 
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Figure 10: The DO 1800 GeV, \rj\ < 0.5 inclusive cross section. Statistical uncer- 
tainties are invisible on this scale. The solid curves represent the ztlo" systematic 
uncertainty band on the data. 

in agreement with the data; in particular, the data above 350 GeV show no indi- 
cation of a discrepancy relative to QCD. The DO collaboration has quantitatively 
compared the data and theory with a test incorporating the uncertainty co- 
variance matrix. The matrix elements are constructed by analyzing the mutual 
correlation of the uncertainties in Fig. |ll| at each pair of Et values. Table Q lists 
values for several JETRAD predictions incorporating various PDFs. Each 
comparison has 24 degrees of freedom. The JETRAD predictions have been fit 
to a smooth function of Ej-- All five predictions describe the \r]\ < 0.5 cross 
section very well (the probabilities for to exceed the listed values are between 
47 and 90%). A similar measurement in the 0.1 < \rj\ < 0.7 interval is also 
well described (probabilities between 24 and 72%). The probabilities calculated 
by comparing the data to EKS predictions for // = (0.25,0.5,1.0) x Ej^'"^^ and 
fi = (0.25, 0.5, 1.0) X E^r^^ are all greater than 57% . Perturbative QCD is in good 
agreement with the data with or without large x enhancements to the PDFs. 
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Figure 11: Contributions to the DO , < 0.5 cross section uncertainty plotted 
by component. 
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Figure 12: Difference between DO data and JETRAD QCD predictions with 
Rsep = 1-3 and fi = EJp"-^ /2 normalized to predictions for the range \r]\ <0.5. 
The bands represent the total experimental uncertainty. 
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Figure 13: Difference between DO data and JETRAD QCD predictions with Rsep 
= 1.3 and fi = EJp""^/! normalized to predictions for the range 0.1< <0.7. 
The bands represent the total experimental uncertainty. 



4.4 Comparison of the CDF and DO Measurements 

The two inclusive jet analyses differ in an important but complementary way. 
The CDF analysis utilizes a Monte-Carlo program carefully tuned to collider 
and test beam data, to correct for detector response. In contrast, DO corrects for 
calorimeter response by direct utilization of collider data. It is worth mention- 
ing that both collaborations invested years of time and effort developing these 
procedures. The CDF technique capitalizes on the excellent tracking capabilities 
of the detector. The tracking is used to directly measure the jet fragmentation 
functions as well as verify test beam calibrations of the calorimeter modules. Ad- 
ditional checks of the Monte Carlo simulation come from comparisons to collider 
jet Et balancing data. The uniformity of the DO calorimetry permits a precise 
collider data-based measurement of jet response and resolution at all energies 
and rapidities. This uniformity permits transfer of the jet energy calibration at 
low to moderate jet energies in the central calorimeter to all pseudorapidity with 
a missing Et technique. The forward calorimetry provides an opportunity for 
direct calibration of high energy jets relative to these central photons. The uni- 
formity, linearity and depth of the detector also assures that resolution functions 
are relatively narrow and Gaussian. 



The top half of Fig. H shows {D - T) /T for the Run IB DO and CDF data sets 
in the 0.1 < \ri\ < 0.7 region relative to a JETRAD calculation using CTEQ4HJ, 



0.5^;^?''=^, 



and 7^sep=l-3. Note that there is outstanding agreement between 



the nominal values for Et < 350 GeV. At higher Et the two results diverge 
but not significantly, given the statistical and systematic errors. This impres- 
sion is fortified by a direct comparison of the magnitude of the D0 and CDF 
uncertainties as shown in the second half of the figure. To quantify the degree of 
agreement, the DO collaboration has carried out a comparison between their 
data and the nominal curve describing the central values of the CDF IB data. 
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Figure 14: Top: comparisons of DO and CDF data to JETRAD in the region 
0.1< \r]\ <0.7. Bottom: the quadrature sums of the D0 and CDF uncertainties. 



The nominal curve is used instead of the data points because each of the mea- 
surements report the cross section at different values of jet Et- Comparison of 
the DO data to the nominal Run IB curve, as though it were theory yields a 
of 41.5 for 24 degrees of freedom. A "statistical-error-only" comparison of the 
DO and CDF data is approximated by calculating the value of the CDF curve at 
the DO Et points and assuming the statistical uncertainty on the CDF and DO 
data are equivalent (the DO statistical errors are multiplied by \/2-) When the 
2.7% relative normalization difference |37] is removed, this statistical-error-only 
is 35.1 for 24 degrees of freedom, a probability of 5.4%. When the systematic 
uncertainties in the covariance matrix are expanded to include both the DO and 
CDF uncertainties and the DO statistical errors are increased by \/2 the equals 
13.1 corresponding to a probability of 96%. 

Considering the complexities of these analyses, the overall agreement is re- 
markable! Given the present flexibility of the PDFs and the mixed agreement 
between the data and theory at the highest x and Q, no clear indication of a 
high Et deviation with QCD can be inferred. By using the high Et jet data, 
PDFs can be derived which describe both data sets. An unambiguous search for 
deviations, however, requires an independent measure of the PDFs. Prospects 
for improvements to the PDFs will be covered in later sections. First we turn our 
attention to inclusive jet measurements at a different beam energy. 



5 RATIO OF JET CROSS SECTIONS AT TWO BEAM EN- 
ERGIES 



Inclusive Jet and Dijet Production at the Tevatron 



25 



5. 1 Inclusive Jet Cross Sections, Scaling, and the Ratio of Dimensionless 
Cross Sections 

An alternative way to test QCD is to compare the inclusive jet cross sections mea- 
sured at widely separated center-of-mass energies. The hypothesis of "scaling" 
predicts that the dimensionless, scaled jet cross section, 

will be independent of as a function of the dimensionless variable xt = 
2E'T/y/s- This can be written in terms of the experimentally measured quan- 
tities as, 

CTrf = {E^/2'K)d^cj/dETdr]. 

QCD predictions depend on the energy scale, or Q^, of the interactions and thus 
suggest that the cross sections should not "scale". The running of the strong 
coupling and the evolution of the PDFs are manifestations of this energy (or 
scale) dependence of the predictions. 

Measurement of the ratio of the scaled jet cross sections from two different 
center-of-mass energies but in the same experimental apparatus, provides a test 
of QCD in which many theoretical and experimental uncertainties cancel. Fig- 
ure |l^ shows the ratio of the scaled cross sections, a^^^ /a^^^^, calculated with 
JETRAD. The top two panels in the figure show a 10% variation in the ratio be- 
low XT = 0.4 (jet Et = 360 GeV at ^=1800 GeV) due to the choice of scale. The 



bottom two panels of Fig. 15 demonstrate that PDF choices produce variations 
below 10%. 

5.2 Jet Production at 630 CeV 

In Dec. 1995 CDF and DO each collected ss600 nb^^ at ^^=630 GeV (beam 
energies of 315 GeV). A primary motivation for this data run was the hope that 
the inclusive jet cross section at this reduced energy and the ratio relative to 
1800 GeV would shed light on the deviations seen with the Run lA 1800 GeV 
cross section. The analysis of the 630 GeV data follows an identical path to the 
analysis described above for the 1800 GeV data. The preliminary measurements 
from CDF p|] and DO pi of the inclusive cross section at ^/s =630 GeV are 



shown in Fig. [L6|. The CDF measurement is for the region 0.1 < |r/| < 0.7 while 
the DO result is for the region \ri\ < 0.5. The figure shows the percent difference 
between the data and the associated theory prediction. The NLO theoretical 
predictions used the MRSA' PDF and = O.Sii^jf*. The two measurements 
are in agreement above 80 GeV, but some discrepancy exists near and below 60 
GeV. The discrepancies are within the 20-30% systematic uncertainties reported 
by D0 and represented by the shaded boxes. With regards to theory, the QCD 
prediction is larger, though not significantly, than the data for Et less than 80 
GeV. Note that in this Et range, results at -v/s=1800 GeV show good agreement 
between data and theoretical predictions. 

5.3 The Scaled Cross Sections 



Figure 17 shows the preliminary ratio of the scaled cross sections Mq] from CDF, 



'^d I'^d ^-long with the previous CDF result [45| for cr^ /c^ from much 
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smaller data samples. The shaded band indicates the systematic uncertainty on 
the 546/1800 ratio; the systematics on the 630/1800 GeV ratio are similar, but 
not yet final. Good agreement between the results is observed. The 546 GeV 
result ruled out "scaling" at the 95% confidence level and a disagreement with 
QCD predictions was observed at low Et- 

Fig. pi shows the preliminary cross section ratios from DO [49| and CDF com- 



pared to a variety of theoretical predictions. Many of the energy scale and lumi- 
nosity uncertainties cancel in the ratio. The uncertainty on the DO ratio, shown 
by the shaded boxes, is about 7%, much less than the 15-30% uncertainty on the 
cross sections, and about a factor of two better than previous ratio measurements. 
The CDF and DO ratios are consistent with each other for xt > 0.1, but some 
difference may exist for xt < 0.1. The discrepancy between the two measured 
ratios is a reflection of the discrepancy in Fig. [l^. The significance of the differ- 
ence must await completion of the CDF systematic uncertainties. The theory is 



roughly 20% higher than the ratio measured by DO. Figure 15 also shows three 
NLO QCD JETRAD predictions for the ratio using /i = O-bE^^"", 7^sep=1.3 and 
different PDFs. In addition, two NLO predictions from the EKS program using 
CTEQ3M and MRSA' and // = 0.5Et are shown. Notice that the variation in 
the predicted ratios is very small. 

The preliminary inclusive jet cross sections at ^/s =630 GeV are not well 
described by NLO QCD calculations. Quantitative results on these comparisons 
await determination of final experimental uncertainties. The ratio of inclusive 
cross sections is also in mild disagreement with the theory. With a larger data 
sample these measurements could place constraints on the high x behavior of the 
PDFs while using relatively low Et jets. An additional run at a similar beam 
energy should be considered for Run 2. We now turn to a discussion of a different 
technique of probing high x behavior: the study of the correlations between the 
leading two jets resulting from a pp collision. 

6 DUET DIFFERENTIAL CROSS SECTIONS AT LARGE RA- 
PIDITY 

6. 1 Introduction 

Measurements of the dijet differential cross sections in different rapidity regions 
can provide additional information and constraints on the QCD predictions. By 
restricting the Et and rapidity of the leading two jets in the events, different 
regions of x and can be probed. This may permit a more direct measure of 
the proton PDFs at Tevatron energies. For instance at LO, an Et = 90 GeV, 
7] = (9=90°) jet must be balanced by a second Et = 90 GeV jet. If the second 
jet is at r/ = 0, both jet energies equal 90 GeV and so xi = X2 = 0.10. However 
if the second jet is more forward, since Et = EsinO and sin^ is smaller, the 
jet energy and its fraction of the initial hadronic momentum must increase to 
maintain Et = 90 GeV. At LO the parton momentum fraction x is related to 
the Et and r] of the two jets by the equations: 

XI = :^(e''i +e''^), X2 = ^(e-^i +e-''^). 
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Figure 15: The difference between alternative predictions and the reference pre- 
diction (/i = O.bEJp""^ , CTEQ3M) of the ratio of inclusive jet cross sections at 
^/s= 630 and 1800 GeV for r/jet < 0.5. The alternative predictions are for the 
choices (a) n = 0.25, 0.75, 1.0 and 2.0 x^f"^, (b) fi = 0.25, 0.5, and 1.0 x^s, 
(c) CTEQ4M, CTEQ4HJ, MRS(Aj, and MRST, and (d) calculations with the 
CTEQ4A series of PDFs (which vary Og) compared with the calculation using 
CTEQ4M. 
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Figure 16: Preliminary D0 and CDF ^/s = 630 GeV cross sections compared to 
NLO QCD predictions. The shaded boxes represent the D0 systematic errors. 
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Figure 17: CDF scaling result from Tevatron runs at \/s= 546 and 1800 GeV in 
1989, compared to the preliminary result from the Run IB runs at 630 GeV and 
1800 GeV. 
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Figure 18: Preliminary D0 and CDF cross section ratios for ^ys = 630 GeV to 
-^s = 1800 TeV compared to NLO QCD predictions. The shaded boxes represent 
the D0 systematic errors. 
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For Et = 90 GeV jets at 77 = and 2, x = 0.06 and 0.42. For multi-jet production 
the calculations of x are generalized to 

xi = ^ V ETie"^' , X2 = ^ V Ene-'^' . 

where "i" runs over all jets in an event. 

Experimentally, the differential dijet cross sections are most conveniently given 

by 

dETdrjidrj2 ' 

where r/i and r]2 are the pseudorapidities of the leading two jets. As with the 
inclusive cross section, the dijet differential cross sections are integrated over a 
range of pseudorapidity. The cross sections are also integrated over a range of 
Et, but here there is even more freedom in the definition of Et'- the leading 
jet Et , the average of the leading two jet E^s, or both leading jet E^s (two 
entries into the cross section, one for the Et of each of the leading two jets). 
The specific choice depends upon the experimental conditions. Although they 
are still preliminary we present the dijet differential cross sections here in order 
to demonstrate the potential and strength of these measurements. Once the 
systematic uncertainties are well understood these measurements will strongly 
constrain the PDFs. 

6.2 The CDF Measurement 

In the CDF Run IB measurement, the two highest Et jets are identified and 
one is required to be in the central (0.1< |?7| <0.7) region. Because the central 
region has the smallest energy scale uncertainty, the central jet is used to mea- 
sure the Et of the event. The other jet, called the "probe" jet, is required to 
have Et >10 GeV and to fall in one of the 77 bins: 0.1< \r]\ <0.7, 0.7< |?7| <1.4, 
1.4< |r/| <2.1, 2.1< |r/| <3.0. There are no restrictions on the presence of addi- 



tional jets. Figure 19 shows the preliminary cross section [^] in the individual 
r] bins as a function of the central jet Et- JETRAD is used for the theoretical 
predictions with scale /i = 0.5 EJp""^ and i?sep=l-3. The data are compared to 
the predictions using three parton distribution functions, CTEQ4HJ, MRST and 



CTEQ4M. The statistical uncertainty is shown on the points. Figure 20 shows 
the percent difference between the data and theory (CTEQ4M) as a function of 
the central jet Et for each rapidity bin. The additional curves represent the per- 
cent difference from the prediction using CTEQ4M for predictions using MRST 
and CTEQ4HJ. The systematic uncertainty on the measurement is evaluated in 
a manner similar to the inclusive jet cross section. The quadrature sum of the 
systematic uncertainties is shown in the box below the points. The high rapid- 
ity bins reach x of roughly 0.6-0.7, while the Et of the jets is in the 100 - 200 
GeV range. The excess over predictions using CTEQ4M observed in the inclusive 
cross section at high Et is seen in all four rapidity bins. This suggests that it 
is not a function of the jet Et, but rather a function of x and is thus related to 
inadequacies in the PDFs. The improved agreement in all bins when CTEQ4HJ 
is used is similarly suggestive. 
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Figure 19: CDF cross sections for central jets with the second jet in different 
rapidity intervals. 
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Figure 20: Percent difference between preliminary CDF data (points) and a 
NLO QCD prediction using the JETRAD program, CTEQ4M, i?sep=l-3 and 



/i=.Ej?"^/2. The percent difference between predictions using CTEQ4M and pre- 
dictions using MRST and CTEQ4HJ are also shown. The four plots represent 
the cross section as a function of the central jet Et when the rapidity of the 
second jet is restricted to different ranges. The curve in the lower box represents 
the quadrature sum of the correlated systematic uncertainties. 
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Table 3: x ranges of SS and OS dijet differential cross sections. 



7] bin 


Topology 




^max 


0.0-0.5 


SS 


0.04 


0.53 


0.0-0.5 


OS 


0.07 


0.42 


0.5-1.0 


SS 


0.03 


0.73 


0.5-1.0 


OS 


0.08 


0.45 


1.0-1.5 


SS 


0.02 


0.81 


1.0-1.5 


OS 


0.10 


0.44 


1.5-2.0 


SS 


0.01 


0.80 


1.5-2.0 


OS 


0.16 


0.46 



6.3 The DO Measurement 

The DO Run IB (92 pb~^) measurement organizes the dijet differential cross 
section according to the rapidity of both leading jets. The rapidities are divided 
into four same-side (SS) bins where 771 ~ rj2 and four opposite-side (OS) bins 
where rji ~ — '?/2- The bins and approximated x ranges sampled (assuming a LO 
process and the observed Et ranges in each bin) are listed in Table |3[ The eight 
cross sections are all plotted versus jet Et- Here each event has two entries, one 
for the Et of each of the leading two jets. As indicated in the table, the low 
rapidity measurements can provide confirmation of previous PDF measurements 
extrapolated from low Q^. On the other hand, the large rapidity, SS events probe 
much larger x values. 

The DO analysis follows the inclusive analysis very closely. Jet and event se- 
lection, energy correction, and resolution unsmearing are all quite similar. An 
additional correction for vertex resolution is also important in the very forward- 
backward rapidity bins. The cross section uncertainties are quite similar to the 
inclusive cross section for the four rapidity bins limited hy rj = 1.0. In the high 

show 
. The 



T] bins the systematic errors are approximately doubled. Figs. |21| and |2^ 
the fractional difference between data and theory for all eight bins [47, [5Q 
theoretical prediction is from JETRAD with CTEQ3M and = 0.5^;^"^. The 
bars on the data represent the statistical uncertainty and the outer symbols the 
total uncertainty. Note there is good agreement over all rapidity for x = 0.01 to 
0.80. Significantly, this PDF includes no collider jet data. Although not shown, 
the agreement with CTEQ4M is also reasonable. 



6.4 Prospects 

The differential cross sections show great promise for constraining the PDFs. 
In all rapidity bins the CDF data appears to prefer CTEQ4HJ over CTEQ4M 
or MRST, while the DO data seems to be good agreement with the predictions 
using CTEQ3M. This apparent disagreement mirrors the situation in the inclusive 
cross sections. However, a firm statement on the agreement or disagreement of 
the two data sets is obscured by the different techniques and the current lack of 
quantitative comparisons between data and theory. 
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Figure 21: DO same-side and opposite-side cross sections compared 
JETRAD with Rsep=l-3, m = E^""^ /2 and CTEQ3M. The two leading 
are restricted to different 77 regions. See the text for details. 
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Figure 22: DO same-side and opposite-side cross sections compared 
JETRAD with Rsep=l.S, n = Ef''''/2 and CTEQ3M. The two leading 
are restricted to different 77 regions. See the text for details. 
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7 DUET MASS AND ANGULAR DISTRIBUTIONS AT 1800 
TEV 

7. 1 Introduction 

The LO cross section for pp jeti + jet2 + X events (where jeti and jet2 are the 
leading two jets) can be completely described in terms of three orthogonal center- 
of-mass variables. These are cos 9*, where 9* is the center-of-mass scattering an- 
gle between the two leading jets, the boost of the dijet system rjboost = {vi + ^2)/2, 
and the dijet mass, Mjj as follows ||5l|| : 

d^<^ ^ ^aim . y fjxi, Q^) f{x2, Q^) 2^ 

dr]f,oostdMjjdcos6* 2s'^ xi X2 ' 

where a1{Q'^) is the strong coupling strength, |mi2p is the hard scattering matrix 
element, xi {X2) is the fraction of the proton (antiproton) momentum carried by 
the parton and f{xi,Q^) is the parton momentum distribution. Typically the 
dijet mass is derived from measured variables such as Et , Vjet, and (j)jet- In the 
case of higher order processes, where more than two jets are produced, the mass 
is calculated using the two highest Et jets in the event and additional jets are 
ignored. 

Integration of the general dijet cross section over boost and production angle 
leads to the dijet mass spectrum. The spectrum is a useful test of QCD sensitive 
to the PDFs. On the other hand integration over mass and boost leads to the dijet 
angular distribution - a marvelous test of the hard scattering matrix elements 
almost totally insensitive to the PDFs. Comparisons of the angular distributions 
and certain mass spectra ratios to theoretical predictions can establish stringent 
limits on the presence of conjectured quark constituents. We turn now to a 
discussion of the these spectra and compositeness limits. 



7.2 The Mass Distributions 

At LO, where only two jets are produced, the dijet invariant mass is given by 

Mjj = s = X1X2S 

where s is the center of mass energy of the interacting partons and s the total 
center ~of-mass energy. Since the dijet mass represents the center-of-mass energy 
of the interaction it directly probes the parton distribution functions. Experi- 
mentally, the dijet mass cross section is given by 

dMjjdr]idr]2 ' 

where t]i and r/2 are the pseudorapidities of the jets. As with the inclusive jet cross 
section the dijet mass cross section is integrated over a range of pseudorapidity. 
For example. Fig. shows NLO QCD dijet predictions for \r]jet\ < 1-0 relative 
to a reference prediction. There are 20-30% variations due to the scale and to 
the PDFs. 

The DO and CDF collaborations have both measured the dijet mass spectrum 
with Run IB data samples. The CDF measurement uses the four-vector defini- 
tion for the dijet mass Mjj = \J {E\ + E2)'^ — {P\ + ^2)'^ where E and P are the 
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Figure 23: The difference between alternative predictions and the reference pre- 
diction {fi = O.bEJp"'^ , CTEQ3M) for the inclusive dijet mass cross section at ^/s= 
1800 GeV for rjjet < 1-0. The alternative predictions are for (a) /j. = 0.25, 0.75, 
1.0 and 2.0 x^^"^, (b) = 0.25, 0.5, and 1.0 x^, (c) CTEQ4M, CTEQ4HJ, 

MRS(Aj, and MRST, and (d) calculations with the CTEQ4A series of PDFs 
(which vary a^) compared with the calculation using CTEQ4M. 



energy and momentum of a jet and allows the rapidity of the jets to extend to 
\r]jet\ < 2.0. To ensure good acceptance, the CDF analysis also imposes a cut on 
|cos^*| < 2/3, where 9* is the scattering angle in the center-of-mass frame. Fig- 
ure 24 shows the CDF preliminary data compared to JETRAD predictions [^]. 
The error bars represent the quadrature sum of the statistical and systematic 
uncertainties. The data and JETRAD predictions using CTEQ4M, = Q.hE^^* 
and Rsep = 1-3 are in good agreement. Figure 25 shows the percent difference 
between the data and the theoretical prediction. The shaded band shows the 
preliminary estimate of the systematic uncertainties. These uncertainties are de- 
rived in a manner similar to the inclusive jet cross section measurement with 
additional contributions for jets outside the central region. The percent differ- 
ence between the default prediction and those using other PDFs are also shown 
and are all consistent with the data given the systematic errors. However, as with 
the inclusive and dijet cross section measurements, CTEQ4HJ seems to provide 
the best agreement with the data in the high Et region. 

DO has performed the dijet mass measurement in two rapidity regions: central, 
\rijet\ < 0.5, and more forward, \rijet\ < 1-0, |Ar/jet| < 1.6 ||5^. The dijet mass 



at DO is defined assuming massless jets: Mjj = 2Ej,Ej,{cosh{Ari) — cos{A(j))), 
where Ar] and A(p are the rapidity and azimuthal separation of the two jets. Fig- 
ure ^ shows the difference between the data and theory divided by the theory. As 
with the DO inclusive jet measurement there is good agreement with theory, well 
within systematic uncertainties. At larger rapidity the data shows a tendency 
to be slightly above the predictions at high Et- Data-theory comparisons 
similar to those described in Section 4 yield large probabilities. As Fig. ^ illus- 
trates there is sensitivity to the PDF choice but unfortunately its significance is 
obscured by the systematic uncertainties. As with the inclusive jet cross section 
these uncertainties are highly correlated. In Fig. 28, the DO \r]\ < 1.0 result is 
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Figure 24: Dijet mass as measured by CDF compared to the NLO prediction 
from JETRAD with CTEQ4M and /j, = £"^"^72 and Rsep=l-3. The error bars 
indicate the quadrature sum of the statistical and systematic uncertainties 




Figure 25: CDF data compared to predictions from the JETRAD program 
for CTEQ4M and /i = EJp"''' /2 and i?sep=l-3 (full circles). Predictions using 
other PDFs and // = £;f"^/2 are also compared to CTEQ4M: MRST (dotted), 
CTEQ4M with n = E^*^^ (dashed), and CTEQ4HJ (sohd). The error bars indi- 
cate the statistical errors. The shaded area represents the combined systematic 
uncertainty. 
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Figure 26: The difference between tlie DO dijet mass data and the prediction 
(JETRAD ) divided by the prediction for \rijet\ < 0.5 and 0.5 < \rjjet\ < 1.0. 
The soUd circles represent the comparison to the calculation using CTEQ3M with 
= 0.5£^r'""^- The shaded region represents the ztla systematic uncertainties. 

compared to the CDF result and to predictions with CTEQ4M. The kinematic 
cuts of the two analyses have a significant overlap; 59% of the CDF sample has 
the two leading jets within \r]\ < 1.0. There is remarkable agreement between 
these data samples over the full Et range. 

7.3 The Dijet Angular Distributions 

The dijet angular distribution measured in the centcr-of-mass is sensitive pri- 
marily to the hard scattering matrix elements. In fact, the distribution is unique 
among high Et measurements in that it is almost independent of the PDFs. 
The shape of the angular distribution is dominated by t-channel exchange and is 
nearly identical for all dominant scattering subprocesses (e.g. gg gg, qg qg, 
and qq — > qq). The dijet angular distribution predicted by LO QCD (two jets 
only) is proportional to the Rutherford cross section 

da 1 

cx 

dcose* sin^{e*/2) ' 

To flatten out the distribution and facilitate the comparison to theory, the vari- 
able transformation % = is used giving 

d^a 

dMjjdxdriboost ' 
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Figure 27: For \rijet\ <1.0, the difference between the DO dijet mass data 
and JETRAD predictions divided by the predictions. The calculations used 
H = O.S^Jr™''"^, Rsep=l-2 and CTEQ4M, CTEQ4HJ, MRS(A'), and MRST. The 
error bars represent the statistical uncertainty on the data. The shaded region 
represents the itlcr systematic uncertainties 
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Figure 28: DO (open circles) and preliminary CDF (full circles) dijet mass 
results compared to predictions from the JETRAD program for CTEQ4M and 
fi = Et/2. 



For measurement of the dijet angular spectra this quantity is integrated over 
pseudorapidity and mass regions and normalized to the total number of events N 
in those regions. The normalization reduces both experimental and theoretical 
uncertainties. Figure 29 shows the jj^^ distributions from the DO collabora- 
tion for four mass ranges |5^. In all cases the jets are limited to regions of 
full acceptance. The NLO predictions from JETRAD with ^ = EJ^'^^ provide 
the best agreement with the shape of the data. The large x acceptance of the 
measurement allows discrimination between LO and NLO predictions. Figure ^ 
shows CDF data compared to QCD predictions (JETRAD ) for different mass 
bins |5^. In this case both LO and NLO QCD are in good agreement with the 
distribution. 



7.4 Compositeness and New Physics Limits 

The dijet mass and angular distributions are sensitive to new physics such as 
quark compositeness. In QCD parton-parton scattering, the dominant exchange 
involves the t-channel. This produces distributions peaked at small center ~of- 
mass scattering angles (near the beam axis in the lab) e.g. large rj and x- 
contrast, the compositeness model predicts a more isotropic angular distribu- 
tion. Thus, relative to QCD predictions, the contributions of composite quarks 
would be most noticeable in the central region, near rj = and x = 1- 

Compositeness signals may be parameterized by a mass scale A which char- 
acterizes the quark-substructure coupling. Limits on A are set assuming that 
A ^ -v/§, such that the dominant force is still QCD. The substructure coupling 
is approximated by a four-Fermi contact interaction giving rise to an effective 
Lagrangian []Tl| . The Lagrangian contains eight terms describing the coupling of 
left and right handed quarks and antiquarks. Currently only the term describ- 
ing the left handed coupling of quarks and anti-quarks has been calculated, and 
this term has an unknown phase. Limits are reported for the case where specific 
quarks or all quarks are composite with either constructive interference (A~) or 
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Figure 29: Dijet angular distribution as measured by DO for different mass ranges 
compared to LO and NLO QCD predictions. 




Figure 30: Dijet angular distribution as measured by CDF for different mass 
ranges compared to LO and NLO QCD predictions. 
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Figure 31: Dijet angular distribution as measured by CDF compared to QCD and 
to QCD plus a term for composite quarks. Limits on compositeness are derived 
from the dijet angular distribution ratio of the number events below x=2.5 to 
the number between 2.5 and 5. 



destructive interference (A"*"). To set compositeness limits CDF and DO both use 
the NLO JETRAD prediction times a "k-factor" from LO QCD + composite- 
ness |54|. NLO calculations with compositeness are not available. Figure pO| 
includes a curve with a compositeness signal added. The additional contribution 
at low X is most pronounced in the highest dijet mass bins. 

Both CDF and DO set compositeness limits based on the ratio of number of 
events at low x to those in the high x region. Figure 31 shows the CDF result 
for the ratio of the number of events below x=2.5 and between x= 2.5 and 
5.0 as a function of the dijet mass, along with curves which correspond to the 
different values of A. Using this ratio, the CDF data excludes at the 95% C.L. 
a contact interaction scale of A^^ < 1.6 TeV and A~^ < 1.4 TeV. For a model 
where all quarks are composite A^ < 1.8 TeV and A" < 1.6 TeV. Figure ^ 
shows the angular distribution as measured by DO for very high dijet masses 
(greater than 635 GeV/c'^) compared to predictions which include composite 
quarks. Compositeness limits in this analysis are derived from the ratio of the 
number of events above and below x=4. The data excludes (at the 95% CL) 
contact interaction scales with n = l.O^f"^: A;;^ < 2.0 TeV, A" < 2.2 TeV, 
A+ < 2.1 TeV, and with ^ = 0.5^f"^: A^^ < 2.2 TeV, A" < 2.4 TeV, A+ < 
2.3 TeV. 

The best limits on compositeness now come from recent dijet mass measure- 
ments from the DO collaboration [^2| which combine the sensitivity of the dijet 
mass distributions with the PDF independence of the angular distributions. In 
much the same way jet production from a composite interaction increases the 
X spectrum at low values of x^ the dijet mass spectrum will increase at central 
rapidities relative to forward rapidities. Thus the ratio of the central and for- 
ward dijet mass spectra will be sensitive to A. In addition, both theoretical and 



experimental uncertainty are reduced in the ratio. Fig. 33 shows the ratio of 
cross-sections for \r]jet\ < 0.5 and 0.5 < \r]jet\ < 1-0 as a function of dijet mass. 
As indicated by the family of curves, the compositeness model predicts changes 
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Figure 32: Dijet angular distribution as measured by DO compared to predictions 
with additional contributions from composite quarks. 



in shape to this ratio at high mass. The spectrum rules out quark compositeness 
at the 95% confidence level for A"*" below 2.7 TeV and below 2.5 TeV. 



8 CONCLUSIONS 

The inclusive jet and dijet cross sections provide fundamental tests of QCD pre- 
dictions at the highest jet and thus are the deepest probes into the structure of 
the proton. With the increased luminosity performance of the Tevatron, measure- 
ments of these cross sections are no longer limited by statistical uncertainties. In 
fact, the systematic uncertainties from the experimental measurements and from 
the theoretical predictions are comparable in size and are significantly larger than 
the statistical uncertainty on all but the highest Et data points. Uncertainties 
in the PDFs dominate the theoretical uncertainty, while uncertainties in the jet 
energy scale dominate the experimental measurements. 

The inclusive jet cross sections from the Tevatron have proved a particularly 
interesting test of QCD. The Run lA inclusive jet measurement showed disagree- 
ment with concurrent pQCD predictions at the highest x and Q. Two subsequent 
measurements, each using five times the data, show mixed agreement between 
the data and theory at high Et- One measurement is consistent with the Run lA 
measurement, and can be described by QCD if the PDFs are suitably modified. 
The second Run IB measurement is well described by QCD with PDFs that either 
do or do not include high Et jet data. Further, within statistical and systematic 
errors, the three measurements are compatible. The preliminary measurements of 
the inclusive jet cross section at v^=630 GeV and the ratio of the cross sections 
show good agreement with previous results and marginal agreement with QCD 
predictions. Derivation of quantitative results is in progress. Unfortunately this 
sample is too statistically limited to provide a constraint on the high Et behavior 
of the cross section at 1800 GeV. 

The apparent excess of events at high Et observed in the inclusive jet cross 
section from Run lA created intense scrutiny of the theoretical predictions. As a 
result, there is now a better understanding of the uncertainty in the theoretical 
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Figure 33: The ratio of DO dijet mass cross sections for \r)jet\ < 0.5 and 0.5 < 
\rijet\ < 1-0 for data (solid circles) and theory (various lines). The error bars 
show the statistical and systematic uncertainties added in quadrature, and the 
crossbar shows the size of the statistical error. 



predictions, particularly for the parton distribution functions. The dijet mass 
distributions and differential cross sections reflect a story similar to the inclusive 
jet measurements. The results are generally consistent with each other and with 
QCD, but with hints that moderately forward regions indicate some high x mod- 
ifications of the PDFs might be appropriate. Incorporation of this information 
into the global fitting procedures used to derive the PDFs should be particularly 
helpful in reducing the uncertainty associated with extrapolation from low energy 
DIS data to the kinematic region covered by the Tevatron jet measurements. 

Until the uncertainty in the PDFs is reduced, or a method of analytically eval- 
uating their uncertainties is available, the best limits on compositeness will come 
from measurements insensitive to the PDFs, and in particular, measurements 
that rely on angular information. These angular distributions show the partonic 
hard scattering to be well described by NLO theoretical calculations. Compar- 
isons of the angular and mass data to jet production models augmented by quark 
compositeness show no preference for deviations from standard QCD predictions. 
In fact, the analyses now show that the compositeness scale must be greater than 
2.5 TeV if it exists at all. 

The high precision Tevatron jet data has fostered a period of great progress in 
QCD. Due to the flexibility of theoretical predictions, pQCD can describe nearly 
all inclusive jet and dijet observations. Limits on quark substructure from Run 1 
have nearly doubled from the previous measurements. However, more stringent 
tests will require improved PDFs and reduced theoretical uncertainties. Looking 
forward, expectations are high that the experimental measurements of jets and 
their properties will continue to improve in Run 2 with a 20-fold increase in data 
sample size, increased beam energy, and reduced systematic uncertainties from 
the upgraded detectors. 
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